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A new numerical method to determine the boundary condition changing (bcc) operators in the 
statistical models is introduced. This method is based on a variant of Schramm-Loewner Evolution 
(SLE), namely SLE(k, p). As a prototype, Abelian Sandpile Model (ASM) with a sink on some 
point on the boundary is considered. Using this method we study the bcc operator corresponding 
to sink point on the boundary. It is numerically shown that the conformal dimension of the this 
operator is nearly 0. The most appropriate candidate for this operator is the logarithmic partner of 
the unity operator / =: 99 :, as it has been conjectured theoretically. 



I. INTRODUCTION 

In determining the operator content of boundary con- 
formal field theory (BCFT), one of the most important 
challenges is to obtain the boundary condition chang- 
ing (bcc) operators corresponding to the changes of the 
boundary conditions. This question has been addressed 
in some minimal models, such as percolatiorP and Ising 
model with spin change on the boundary^ etc. To de- 
termine the conformal weight of these bcc operators, one 
has to determine the correlation functions of the model 
and compare them with the exact results, or use some ev- 
idences e.g. operator product expansion (OPE) of these 



operators^. This needs almost complete knowledge of 
the operator content of the model in hand which for the 
less-known models is not applicable. Therefore one needs 
more direct method to study such operators. An impor- 
tant recent discovery, concerning this difficulty, was the 
discovery of the bcc operator corresponding to the change 
of boundary conditions of the points at which the trace 
of a critical interface starts. These interfaces are some 
non-intersecting macroscopic objects which separate the 
different phases of the model and are postulated to be 
described by Schramm-Loewner Evolution (SLE) theory 
whose operator content is almost well understoocP. 

Schramm-Loewner evolution, as the new interpreta- 
tion of statistical models, has attracted much attention 
in last years. According to this theory one can classify 
2D statistical models via a growth processes by focusing 
on their geometrical objects, such as interfaces. The es- 
sential building block of SLE is the conformal symmetry 
of probability measures of the model in hand. Therefore 
it seems natural to exist a connection between SLE and 
conformal field theory (CFT) in which, contrary to SLE, 
one deals with the local fields. This relation shows it- 
self in a simple relation between the central charge c in 
CFT and the diffusivity parameter k in SLEP. The pa- 
rameter k is a quantity which represents the universal- 
ity class in which the statistical models belong to. Due 
to its capability in reflecting essential properties of the 
models, SLE has been referred and employed theoreti- 
cally and numerically in various statistical models such 
as turbulence^, ASM avalanche frontier^, iso height lines 
of KP#>, WOS^, Ising modeP, etc. to study statisti- 
cal properties of them. Equally much theoretical works 



have been carried out on this theory, such as left passage 
probability and fractal dimension of the curves^, crossing 
probability, etc. 

In this paper we introduce a new method to determine 
numerically the conformal weight of the bcc operator us- 
ing a generalization of SLE. This generalization, namely 
SLE(k, p) governs the law of the curves whose growth me- 
dia contains, in addition to the starting and ending points 
of the curve, some other preferred points. In this theory 
p has to do with the conformal weight of the correspond- 
ing bcc operator and has the capability to yield some 
informations about this operator. We test this method 
on the Abelian sandpile model (ASM) to determine the 
conformal weight of the sink points (at which the grains 
of the sand pile dissipate) on the boundary of the domain 
in which the model is defined. Although there are some 
theoretical evidences that in this model, the bcc opera- 
tor is the logarithmic partner of the identity operator i. 
c. I =: 98 : (z), since ASM is somehow problematic, we 
believe that it is interesting to study and investigate this 
operator directly using our method. 

The next section is devoted to introduction to SLE and 



its generalization SLE(k, p). In section III we introduce 



in sum the ASM and define our main concern in this 



paper and in section IV we present the numerical results 



of the application of our method mentioned above. 



II. SLE 



In two dimensions, the phase boundaries of the statis- 
tical models form some non-intersecting random curves 
which directly reflect the status of the system in ques- 
tion. This theory is applicable to the curves which are 
supposed to have two properties: conformal invariance 
and the domain Markov property. According to this the- 
ory, each conformal invariant statistical model fall into 
a one-parameter universality class which is represented 
by k. We offer a very brief introd uctio n below; for good 
introductory reviews see reference d 6 ! 15 !. 
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A. Chordal SLE 

Let us denote the upper half-plane by H and 7 as the 
SLE trace. Then parametrize the curve with "time" 
t and name the curve up to this time as 7*. Define 
t z as the time at which this curve meets the point z 
in the upper half plane. The hull of this growth pro- 
cess is defined as K t — {z G H : t 2 < t} so that the 
complement H t := H\K t is simply-connected. Accord- 
ing to Riemann's mapping theorem there is always a 
unique conformal map gt{z) which maps H t — > H with 
gt(z) — z + — + 0(\) as z — > 00 known as hydrody- 
namical normalization. SLE K is then defined via such 
conformal maps which are solutions of Loewner's equa- 
tion: 



d t gt(z) = 



(l) 



with the initial condition gt(z) — z. In this equation the 
driving function £ t is proportional to a one dimensional 
standard Brownian motion i.e. £ t = ^J~k~B t . For fixed z, 
gt(z) is well-defined up to time t z for which <?t(z) = £4. 
One can retrieve the SLE trace by 74 = lim e ^o 5t _1 (?t+* e )- 
There are phases for these curves, for < k < 4 the trace 
is non-self-intersecting and it does not hit the real axis; 
in this case the hull and the trace are identical: K t — "ft- 
This is called "dilute phase". For 4 < n < 8, the trace 
touches itself and the real axis so that a typical point is 
surely swallowed as t — > 00 and K t ^ 7 t . This phase 
is called "dense phase". And for k > 8 the curve is 
space-filling. The frontier of K t , i.e. the boundary of H t 
minus any portions of the real axis for 4 < n < 8 is a 
simple curve which is locally a SLE S curve with k = — 
(2 < k < 4), i.e. it is in the dilute phased This relation 
relates two phases to each other. 



hydrodynamically normalized, then the stochastic equa- 
tion of gt is the same as Eq. ([I]). In fact this procedure 
leaves the Eq. ([I]) unchanged but leads the driving func- 
tion to have a drift terrrPIl 



k — 6 



6 - gtixoo) 



dt 



(2) 



Thus for the critical curves from boundary to boundary, 
the corresponding driving function acquires a drift term. 
This generalization of SLE can be generalized further to 
have mul tiple preferred real axis points. For review see 
reference d 6 * 15 ^. Eq[2] is especial case of the more general 
theory i.e. SLE(k, p = k — 6). 

In the SLE(k,p), the parameter k identifies the local 
properties of the model in hand and corresponds directly 
to the central charge of the corresponding conformal field 
theory and the parameter p has to do with the boundary 
conditions (be) imposed i.e. some information of the bec 
operator have been coded in this parameter. The actual 
behavior depends on the concrete values of k and p. The 
example is the dipolar SLE(k) in which p = (k — 6)/sH 
The other example is the situation in which except the 
origin and the ending point (say infinity), the boundary 
condition changes at the other point x$. The stochastic 
equation governing such curves is the same as formula 
([T]) but the driving function has a different form: 



P 



&-gt(xo) 



dt 



(3) 



The generalization of Eq [3] is direct for more preferred 
points on the real axis (xq,xi, ...,x n ' 



120. 



i=0 



&-gt(xi) 



dt 



(4) 



B. SLE(/t,p) 

SLE theory describes the critical random curves which, 
except the origin (at which the curve starts to grow) and 
the ending point (which is infinity in the chordal case), 
there is no other preferred points on the real axis. But 
in some situations the existence of some preferred points 
on the real axis does affect the growth of the curve. As 
an example let us consider a curve that starts from the 
origin and ends on a point on the real-axis (xoo). In 
this situation, we will have two preferred points on the 
real axis. Using the conformal map 4> = x 00 z/(x 00 — z), 
one can send the end point of the curve to the infinity. 
In this respect, the function h t = (jjogt ^ 1 describes 
chordal SLE. It is easy to show that the equation govern- 
ing h t is d t g t = 2/({(j)'(g t )((j)(gt) - &)})■ But it is explicit 
that this function is not hydrodynamically normalized. 
It has been showrf^l that if one uses another mapping 
g t = v t oh t ou~ x where u — (j)" 1 and v t is a linear frac- 
tional transformation that make the corresponding map 



III. INTRODUCTION TO ASM AND THE 
PROBLEM DEFINITION 

The simple model of Bak, Tang and WiesenfeloP^ 
of Self-Organized Criticality (SOC) phenomena has at- 
tracted much attention due to its rich structure and com- 
plex behaviors. As an example of this phenomena, they 
introduced the sandpile models in which without tuning 
any parameter, the system show critical (such as power 
law) behaviors. The abelian structure of this model was 
first discovered by D. Dhar and named as Abelian Sand- 
pile Model (ASM)P9 ASM has various and interest- 
ing features and many different analytical and numerical 
works have been done on this model. For example differ- 
ent height and cluster probabilities 25 , its connection with 
spanning treea^, ghost models^, q-state Potts modeP^, 
etc. For a good review see referenced 

Before going to a specific geometry, let us consider the 
ASM on a two-dimensional square lattice Lx L. To each 
site i, a height variable hi is assigned taking its values 




FIG. 1: The schematic picture of the Abelian sandpile model 
on the cylinder in which the grey region has toppled. 



from the set 1,2,3,4 the number of sand grains on this 
site. Define the dynamics of this model as follows; in each 
step, a grain is added to a random site i i.e. hi — > hi + 1; 
if the resulting height becomes more than 4, the site top- 
ples and loses 4 sand grains, each of which is transfered 
to one of the four neighbors of the original site. As a 
result, the neighboring sites may become unstable and 
topple and a chain of topplings may happen in the sys- 
tem. At the boundary sites, the toppling causes one or 
two sand grains leave the system. This process continues 
until the system reaches a stable configuration. Now an- 
other random site is selected and the sand is released on 
this site and the process continues. After a finite number 
of steps, the ASM reaches a well-defined distribution of 
states in which some configurations do not occur (tran- 
sient configurations) and other configurations occur with 
equal probability (recurrent configurations). For details 
see reference^. This model is related to Potts model 
with q — > and CFT with central charge c = —2. For 
a lattice with d neighboring sites, the toppling occurs 
when hi > d, then the original site will lose d grains and 
the height of each of its neighbors will increase by 1. It 
has been shown that for such a system, the avalanche 
frontier is a random loop with the fractal dimension | 
and the same statistical properties as the loop erased 
random walk (LERW)P^. To determine the diffusivity 
parameter k we should find some random curves which 
goes to a point at infinity whereas the mentioned fron- 
tier forms a closed loop. To make such curves, one can 
cut the loops horizontally and send the end point of the 
curve to the infinity. This method causes some large nu- 
merical errors. It has been proposed irP^l that one can 
skip the mapping to infinity and interpret the resulting 
curve as a curve which start and end on the real axis. In 
this situation the formalism of SLE(k,p) can be applied 
to find the diffusivity k and p parameters of the model. 
To this end the uniformizing map firstly should become 
disceretized G tn ,£ n = Gst n _ 1 ^ n _ 1 oGgt n _^ n _ 3 o...oGs to> ( , 
(St n — t n — t n -\) in which 



(5) 



and then assuming the driving function partially con- 
stant (in each time interval), so the Eq [2] can become 
discretized as follows: 



/Ji5B n 



St,. 



(6) 
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FIG. 2: A sample curve which is the cut exterior boundary 
of an avalanche. 



Rewriting this equation in the form: 

t _ r pSU i 



-Cm 



(7) 



one can find the corresponding k and p by demanding 
that the right hand side of Eq [7] be a one dimensional 
Brownian motion. It has been shown that this method 
results in more precise and reliable determination of these 
parameters^. 

Now consider one additional preferred point on the real 
axis sitting in Xq. The random curves in this media are 
described by SLE(k, k — 6, p) in which the parameter p is 
related to the bcc operator in xq. Then the Eq [4] yields: 



FndBt 



6 



-dt 



tt-gt(xoo) it-9t{x ) 



dt. 



(8) 



There is a simple relation between the conformal weight 
of bcc operator at xq, h p and p i.e. h p = p( - p + ±~ K ^ ■ All 
the arguments mentioned above are applicable to this 
case i.e. one can determine the unknown parameters of 
the model by fitting the following equation by one dimen- 
sional Brownian motion: 



[ (Pc)Stj 

=i L ei-G ti (x 



1 — V" r P St i i 

)J l^i=l[^-G t .(x )i 



= B n . (9) 



The best value of p is such a parameter which results 
in a nice one-dimensional Brownian motion in right 
hand of Eq (9). Since the numerical calculation of these 
three parameters togather is difficult, the best method 
to obtain this quantity is to obtain best values of k and p 
for the case there is no boundary condition change, and 
then repeat the calculations for the case with boundary 
condition change, setting k and p c fixed. 
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FIG. 3: The fractal dimension of the curves in two cases; with 
and without sink point. 



FIG. 4: The graph of (B t 2 ) defined in Eq[9] versus t for the 
ASM without sink point (ft = 2.0 and p c — —3.2) and the case 
involving sink point with the same k and p c and p — —0.1. 



A. Fractal dimension 



IV. NUMERICAL METHODS AND RESULTS 



Consider ASM on a cylinder with length and circum- 
ference L in presence of a sink point in which the grains 
dissipate (the point xq). Fig [I] shows schematically such 
a model in which the grey area is the set of toppled sites. 
Therefore in one direction the boundary conditions are 
periodic in contrast to the other direction in which the 
sand grains may leave the system i.e. with open bound- 
ary condition. This geometry is compatible with our pur- 
pose i.e. we have real boundary at which the sink point 
has been located. Now we consider the upper portion 
of the loop (the frontier of the avalanche) in such a way 
that the resulting curve goes from origin to the point Xoo 
on the boundary in presence of the the sink point at xq. 
If L be much larger than the linear size of the curves, we 
can approximately ignore the finite size effect and apply 
the chordal SLE formalism. We also rescale the curves in 
such a way that xq — £o = 1. In figure[2]we have presented 
a sample curve which is the cut exterior boundary of an 
avalanche. IrP by sending the end point of such a curve 
to infinity using making chordal SLE, it was shown that 
the amount of k for this problem is nearly 2 (without 
dissipation in sink). In this paper we do not apply such 
a infinity map due to the reasons to be mentioned in the 
following subsections. 

The simulation has been done over 3 x 10 4 samples at 
which the minimum size of the curves is 1000 in lat- 
tice units. The CPU time for this simulation was nearly 
3 x 10 6 s for each case^. 



To avoid the possible errors due to existence of couple 
of equations in the SLE(k, p c , p) formalism, we have made 
two parallel simulations with the same conditions, one for 
the case in which there is no sink point, and another for 
the case involving the sink point. One can calculate the 
amount of p by comparing the resulting Brownian motion 
of the two simulations. First of all, we should determine 
the best values for k and p c for the first case and then 
fixing these quantities, we can find the best value of p 
for the second case for which the two simulations best 
coincide. To this end we have calculated the fractal di- 
mension of the curves in each case. We use the fact that 
(showing the linear size of the curve by d and the length 
of the curve by I) I ~ d D f where Df is the fractal di- 
mension of the curve. The Fig [3] in which the log-log 
plot of I and d is sketched, demonstrates that the frac- 
tal dimension of the curves in two cases coincides and is 
nearly |. The main result in this subsection is that the 
amount of k and p c in the two cases is equal since the 
fractal dimension is only related to the amount of k i.e. 
Df = 1 + | and p c only depends on «P. 



B. Calculation of p 



To use the result of the subsection IV A[ we should ob- 
tain k and p c for the case without sink point, and then 
use these parameters for calculationg p for the case in- 
volving sink point. The first simulation has been done 
irP^ in which it was shown that the best values are 
k = 1.95 ± 0.07 and p c = 3.5 ± 0.5. Our simulation 
for our case (which is on cylinder and the curves are 
rescaled) coincides with this report, i.e. we have obtained 
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k = 2.0 ± 0.1 and p c — 3.2 ± 0.5. The most important 
part of our results is the second part of simulation, i.e. 
determination of p. Fig [4] shows the variance of Brown- 
ian motion for the two mentioned cases. We see that the 
two curves coincide by attributing to p the value —0.1. 
By using Maximum Likelihood Estimation (MLEpHl we 
have obtained p — —0.1 ± 0.2 with the probability 0.8. 
From this result one can calculate the conformal weight 
of the bcc operator corresponding to sink; h p = O.OIq'q?. 

Our result is in agreement with the prediction of field 
theoretical point of view of ASM. It has been conjectured 
that the sink point, which can be interpreted as a change 
of boundary conditions from closed to open and imme- 
diately from open to closed, is equivalent to putting an 
operator resulting from the OPE of two twist operators p 
with the conformal weig ht -|. According tcP 1 , the OPE 
of p operator is the direct sum of identity operator and 
its logarithmic partner and the second one corresponds 
to the sink point which is our case. So the result of this 



OPE is nothing but the logarithmic partner of identity 
operator I —: 99 : (z) with the conformal weight h p = 0. 



V. CONCLUSION 

In this paper we have proposed a framework to ob- 
tain numerically the conformal weight of the boundary 
condition changing (bcc) operator in a critical statistical 
model. As an example, Abelian sandpile model (ASM) 
has been considered. Specially we have considered a sink 
point on the boundary in which the grains dissipate and 
have calculated the conformal weight corresponding to 
the bcc operator to be h p = O.O^q'qj which is in agree- 
ment with the field theoretical result h p = correspond- 
ing to the conformal weight of logarithmic partner of the 
identity operator I =: 99 : (z). 
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